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1. A method for predicting a plurality of surface multiples for a plurality of traces in 
a record of seismic data, comprising: 

5 providing a plurality of target traces at a nominal offset and a nominal azimuth; 

selecting a plurality of pairs of input traces, wherein the midpoints of the input 
traces in each pair are separated by half the nominal offset and the azimuth of a line 
connecting the midpoints of the input traces in each pair is equal to the nominal 
azimuth; 

10 convolving the selected pairs of input traces to generate a plurality of 

convolutions; and 

applying a three dimensional operator to the convolutions. 

2. The method of claim 1, further comprising preconditioning the input traces to 
simulate zero offset traces. 

15 3. The method of claim 1, further comprising preconditioning the input traces to 
simulate zero offset traces by a moveout correction. 

4. The method of claim 1, wherein the target traces provide the locations at which 
the surface multiples are predicted. 

20 

5. The method of claim 1 , wherein the nominal offset is a central offset value within 
a range of offsets. 

6. The method of claim 1, wherein the nominal azimuth is a central azimuth value 
25 within a range of azimuths. 

7. The method of claim 1 , wherein selecting the pairs of input traces comprises 
selecting the input traces at half the nominal offset. 

30 B. The method of claim 1 , further comprising interpolating and regularizing the input 
traces. 
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9. The method of claim 1 , further comprising locating each convolution at the 
midpoint of the line connecting the midpoints of the input traces in each pair. 

1 0. The method of claim 1 , wherein applying the three dimensional operator to the 

5 convolutions comprises locating the result of the three dimensional operator application 
at the midpoints of the target traces. 

1 1 . The method of claim 10, further comprising correcting the result to an actual 
offset and an actual azimuth. 

10 

12. The method of claim 1 , wherein the three dimensional operator is a three 
dimensional demigration operator having an offset equal to half the nominal offset and 
a velocity equal to one of half the water velocity or half of a multiple velocity function, 

15 13. A method for predicting a plurality of surface multiples for a plurality of traces in 
a record of seismic data, comprising: 

dividing a plurality of target traces into one or more groups according to offsets; 

dividing each group into one or more subgroups according to azimuths; 

selecting a first subgroup having a first nominal offset and a first nominal 
20 azimuth; 

selecting a plurality of pairs of input traces, wherein the midpoints of the input 
traces in each pair are separated by half the first nominal offset and the azimuth of a 
line connecting the midpoints of the input traces in each pair is equal to the first nominal 
azimuth; 

25 convolving the selected pairs of input traces to generate a plurality of 

convolutions; and 

applying a three dimensional operator to the convolutions. 

14, The method of claim 13, wherein the target traces within each group are at a 
nominal offset. 

30 15. The method of claim 14, wherein the nominal offset is a central offset value 
within a range of offsets. 
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1 6. The method of claim 13, wherein the target traces within each subgroup are at a 
nominal offset and a nominal azimuth. 



17. The method of claim 15, wherein the nominal azimuth is a central azimuth value 
5 within a range of azimuths. 

18. The method of claim 13, further comprising preconditioning the input traces to 
simulate zero offset traces. 

19. The method of claim 13, further comprising locating each convolution at the 
10 midpoint of the line connecting the midpoints of the input traces in each pair, 

20. The method of claim 13, wherein applying the three dimensional operator to the 
convolutions comprises locating the result of the three dimensional operator application 
at the midpoints of the target traces. 

15 

21 . The method of claim 1 3, wherein the three dimensional operator is a three 
dimensional demigration operator having an offset equal to half the first nominal offset 
and a velocity equal to one of half the water velocity or half of a multiple velocity 
function. 

20 

22. A method for predicting a plurality of surface multiples for a plurality of traces in 
a record of seismic data, comprising: 

providing a plurality of target traces at a nominal offset; 

predicting a plurality of two dimensional surface multiples for a plurality of input 
25 subsurface lines; and 

applying a two dimensional operator to the predicted two dimensional surface 
multiples. 

23. The method of claim 22, further comprising, prior to predicting the two 
30 dimensional surface multiples, preconditioning the input traces into the input subsurface 

lines, 
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24. The method of claim 22, further comprising sorting the predicted two 
dimensional surface multiples into one or more planes according to offsets, wherein 
each plane has an associated nominal offset. 

25. The method of claim 24, further comprising selecting a plane of predicted two 
dimensional surface multiples having a first nominal offset. 

26. The method of claim 22, wherein applying the two dimensional operator 
comprises locating the result of the two dimensional operator application at the 
midpoints of the target traces. 

27. The method of claim 22, wherein the two dimensional operator is a two 
dimensional demigration operator with a velocity equal to one of half the water velocity 
or half of a multiple velocity function. 
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